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ABSTRACT 

The recently revised cosmological constraints from the Five- Year WMAP data ameliorate previous 
tension between cosmological constraints from the microwave background and from cluster abun- 
dances. We demonstrate that the revised estimates of cosmological parameters are in excellent agree- 
ment with the mass function of X-ray clusters in the Sloan Digital Sky Survey. Velocity segregation 
between galaxies and the underlying dark matter could cause virial mass estimates to be biased, 
causing the mass scale of the mass function to be offset from the true value. Modest velocity segrega- 
tion ((Tg X y/aDM = ^-^-o'ot) is sufficient to match the mass function to the Five- Year WMAP results. 
When the new WMAP results are combined with constraints from supernovae and baryon acoustic 
oscillations, there is no need for velocity segregation (cr g:cy /<7DM=1.05±0.05). This result agrees with 
expectations for velocity segregation from state-of-the-art numerical simulations of clusters. Together 
with the improved agreement between the new WMAP results and recent cosmic shear measurements, 
this result demonstrates that the amplitude of large-scale structure in the nearby universe matches 
that predicted from the structure seen in the microwave background. The new constraint we place on 
velocity segregation in clusters indicates that virial mass estimates for clusters are reasonably accurate. 
This result suggests that future cluster surveys will be able to probe both cosmological parameters 
and fundamental cluster physics. 

Subject headings: galaxies: clusters — galaxies: kinematics and dynamics — cosmology: observations 



1. INTRODUCTION 

The increased statistical power and sophistication of 
different cosmological probes enables new tests of long- 
standing astrophysical problems. We compare recent 
cosmological constraints from observations of the mi- 
crowave background with those obtained from cluster 
abundance estimates. This analysis effectively compares 
the amplitude of large-scale structure at z«1100 with the 
large-scale structure in the local universe. 

Clusters of galaxies are the most massive gravita- 
tionally relaxed systems in the universe, so the ob- 
served cluster mass function is a sensitive probe of cos- 
mological parameters. Galaxy cluster abundances are 
most sensitive to the matter density of the universe 
TO and erg, the rms fluctuations in spheres of radius 
8ft. _1 Mpc and the normalization of the linear powe r 
spectrum (|Henrv fc Arnaudlll99H TBahcall fc Cenlll993fl . 
The estimated values of these two parameters underwent 
significant revisions between the One- Year and Three- 
Year WMAP resul ts (hereafter WMAP1 and WMAP3; 
ISpergel et al.| [2007'l . With the recent release of the Five- 
Year WMAP results (hereafter WMAP5), the estimates 
of these two parameters shifted again, although this 
latest shift is smalle r than the statistical uncertainties 
(|Dunklev et al.ll2008l see their Table 2 and Figure 6). 

The shift in these parameters between WMAP1 and 
WMAP3 significantly alters the expected abundance 
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of massive clusters. iReiprichl (|2006l ) noted that this 
revision agrees well with th e X-ray mass function of 
iReiprich fc Bohringerl (|2002l ). However, other studies 
suggest that h ydrostatic X-ray mass estima tes such as 
those used in IReiprich fc Bohringerl (|2002f ) underesti- 
mate true cluster masses b ecause either the gas temper- 
atures are underestimated (Rasia et al. 2005T) or because 
turbulent flows and kinetic pressure in the intraclus- 
ter medium (ICM) ar e a significant ene rgy component 
(jVikhlinin et al.l l2006t iNagai et al.ll2007h . The (poorly 
constrained) amount of energy in non-thermal pressure 
has a dramatic impact on estimates of eg i nferred from 
X-ray observations (e.g., see Figure 17 of iMantz et al.1 

We recently estimated the cluster mass function using 
ROSAT X-ray cluster surveys to determine the selection 
functio n and redshifts from th e Sloan Digital Sky Survey 
(SPSS. IStoughton et al.l l2002ft to estimate virial masses 
(|Rines et al.ll2007D . The cluster masses were computed as 
part of the Cluster Infa ll Regions in SDSS project (CIRS, 
iRines fc Diaferid 12006) . Our estimate of cluster abun- 
dance exceeded that expected for the best-fit WMAP3 
parameters. We showed that this difference could be at- 
tributed to velocity segregation: if cluster galaxies have 
a significantly higher velocity dispersion than the under- 
lying dark matter, then dynamical mass estimates based 
on the virial theorem overestimate cluster masses and 
thus overestimate the abundance of clusters at a given 
mass threshold. We argued that significant velocity seg- 
regation (ag X y/aDM~l-2&±0-06) is needed to produce 
agreement with th e cosmological parameters of WMAP3 
(|Rmes et al.ll2007t) . 

Numerical simulations including models of galaxy for- 
mation typically predict little velocity segregation (also 
called "ve locity bias " ) in clu sters (|Ka uffman n et all 
fT999alrbl : iGao et al.l l200l iDiemand et all 120041 : 



Faltenbacher et al.ll2005t iFaltenbacher k piemandll2006l: 



Biviano et all 120061: iBenatov et all 120061 : lEvrard et al l 
20081 ). lEvrard et all (|2008l) showed that velocity disper- 



sions of dark matter halos provide robust and accurate 
mass estimates over many orders of magnitude in halo 
mass; that is, dark matter particles obey virial scaling re- 
lations. Further, they confirmed that the C1RS velocity 
dispersion function requires significant velocity segrega- 
tion in clusters to match the cosmological parameters 
estim ated from a joint ana lysis of WMAP3 and SDSS 
data l|Tegmark et alj [2006) . In particular, they show 
that the parameter S$ = as(Q m /0-3) ' 35 has conflicting 
estimates of S 8 =0.69 from WMAP3+SDSS versus 
5*8~0.9 from the CIRS velocity dispersion function. 
lEvrard et all (pool also show that the WMAP3+SDSS 
cosmology requires excess specific energy in the ICM, 
again conflicting with expectations from numerical 
simulations. They propose that a plausible intermediate 
value of 58=0.8 would reconcile the estimates without 
requiring excess specific energy in galaxies or the ICM. 

In this Letter, we show that the revised cosmologi- 
cal constraints from WMAP5 and other methods signif- 
icantly alleviate the need for velocity segregation and 
excess ICM specific energy in galaxy clusters. We re- 
view the challenges of accurately determining the clus- 
ter mass function in §2 and compare the WMAP5 re- 
sults with the cluster mass function in §3. We discuss 
the implications of these comparisons in §4. We assume 
H = 100k km s^Mpc" 1 , and a flat ACDM cosmology 
(£1a = 1 — ^m) throughout. Where not stated explicitly, 
we assume VL m = 0.3 and h=0.7 for initial calculations. 

2. DETERMINING THE CLUSTER MASS FUNCTION 

There are many challenges to accurately estimating the 
abundance of galaxy clusters. The first challenge is accu- 
rately measuring the survey volume. In particular, one 
must know the maximum volume V max in which a clus- 
ter with certain properties could be detected within a 
given survey. Most clusters exhibit X-ray emission from 
the ICM. As a result, X-ray surveys can be used to de- 
fine flux-limited samples where V max can be computed 
directly. Optical surveys are excellent for detecting clus- 
ters, but they are more sensitive to projection effects 
than X-ray surveys. Our approach was therefore to com- 
bine the virtues of a well-defined selection function avail- 
able from large-area X-ray surv eys with the large redshift 
samples available from SDSS (|Rines et al.l I2007T ). This 
approach has the advantage that cluster selection and 
cluster mass estimates are decoupled. 

Another significant challenge is sampling a large and 
representative volume. Clusters are extremely rare ob- 
jects, so a large survey volume is required to fairly 
sample the cluster mass function. The cluster sam- 
ple w e use is the CIRS sample of iRines k Diaferiol 
(120061) with some minor modifications described in 



IRines et al.1 (|2007t ). The CIRS clusters are selected from 
X-ray cluster cat alogs constructed fr o m ROSAT All- 
Sky Survey data (|Ebeling et al.l 120001 : iBohringer et all 
120001 120011 ). We searched for all z<0.1 X-ray clusters 
within the 4783 square degrees of spectroscopic data 
available in the Fourth Data R elease (DR4) of SDSS 
(jAdelman-McCarthy et ai1l2006! ). The volume contained 
in the CIRS mass function is ~10 7 fe 3 Mpc~ 3 , the larges t 
ever probed by virial mass estimates (|Rines et al.ll2007t ). 



The CIRS clusters are an unbiased sample: the selec- 
tion of the sample is based purely on X-ray flux and the 
footprint of the SDSS DR4 spectrosco pic survey. We con- 
firmed this claim with a V/V max test rtSchmidtlfl968f) : we 
find < V/V m ax >=0.518±0.035 compared to an expected 
value of 0.5 for a complete, uniform sample (jRines et al.l 
l2007f h The CIRS sample is essentially volume- limited for 
clusters with Ljf>3xl0 43 /i~ 2 erg s _1 and z<0.1. 

The greatest obstacle to measuring the cluster mass 
function is obtaining sufficiently accurate mass estimates. 
It is obscrvationally challenging to obtain detailed mass 
estimates for large samples of clusters. One common ap- 
proach to computing the cluster mass function is there- 
fore to use relatively simple obs ervables such as X- ray 
temperature or luminosity (e.g., IMantz et al.l I2008D or 
optical luminosit y or richness (e.g.. iBahcall et al.l 120031 : 
IRozo et al.l l2008f as proxies for detailed mass estimates. 
Because these mass proxies require significantly less data 
than detailed mass estimates, they can be applied to 
large samples of survey-quality data and thus probe 
larger volumes than CIRS. However, an important as- 
sumption in this method is that residuals from the scal- 
ing relation must be well-behaved (the residuals are usu- 
ally assumed to f ollow a log- normal distribution, e.g., 
IMantz et ai]|2008D . 

For the CIRS mass function, we use detailed dynamical 
mass estimates for the individual clusters to compute the 
mass function. The dynamics of cluster ga laxies provided 
the fi rst evidence for dark matter when IZwickvl (| 19331 . 
Il937| ) applied the virial theorem to the Coma cluster. 
Since then, the use of virial mass estimators has been 
refined and this method has proved to be a powerful too l 
for measuring cluster m asses (e.g.. iBiviano et aT1l2006h . 
IRines k Diaferiol (|2006l ) show the infall patterns of the 
72 CIRS clusters and compute the mass profiles from 
both the caustic technique and the virial theorem. With 
the dense redshift samples available from SDSS, cluster 
members can be cleanly separated from foreground and 
background galaxies, and the statistical uncertainties in 
the virial mass estimates are relatively small. 

Compared to X-ray studies, virial mass estimates are 
sensitive to larger scales (r2oo rather than rsoo, where ta 
is the radius within which the enclosed density is A times 
the critical density). This difference allows comparisons 
with theoreti cal mass func tions with significantly less ex- 
trapolation (| White! [2002D . Also, virial masses can be 
estimated for poor clusters and rich groups, whereas X- 
ray mass estimates of these systems are complicated by 
possible energy in put from supernovae and AGN (e.g., 
iLoewensteinl 12000). Probing these smaller masses en- 
ables a direct constraint on fluctuations on the scale 
8/i -1 Mpc, rather than the ^14/t~ 1 Mpc scal e probed by 
~10 15 /i" 1 M Q clusters (jPierpaoli et al.ll2001l ). 

We used the virial mass function to obtain cosmo- 
logical constraints; these constraints can be parameter- 
ized as ^(^/O.S ) 5 = 0.81^;^ ( 6 8% co nfidence level; 
IRines et al.l I2007D . lEvrard et all (|2008D showed that 
the CIRS velocity dispersion function yields similar re- 
sults to the CIRS mass function. Using stacked clusters 
from the maxBCG sample to determine scaling relations, 
iBecker et al.l ()2007l ) found a velocity dispersion function 
in good agreement with that found by CIRS. These com- 
parisons support our claim that the CIRS mass function 
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Fig. 1. — The mass function of the CIRS sample. Green and 
blue points and errorbars are computed using virial masses and 
caustic masses respectively (the errorbars show 68% uncertain- 
ties). The dash-dotted lines show the mass functions computed 
using the cosmological parameters from the WMA P1 results (up- 
per) a nd WMAP3 results (lower) using the results of ljenkins et al.l 
(2001). The red line shows the best-fit mass function for the CIRS 
virial mass function. The orange (solid) and magenta (dashed) 
lines show the expected mass functions from the WMAP5 and 
WMAP5+SN+BAO parameters respectively. 

is accurate. 

3. COMPARING TO WMAP5 AND WMAP5+SN+BAO 



Figure [T] shows the mass functions for t he best-fit cos- 
molo gical parameter s from WM A P1 (jSpergel et"aTI 
[2001. WMA P 3 flSpergel et alJ |2Q0J), WMAP5 
([Dunklev et "all 12008^ and WMAP5+SN+BAO 
(|Komatsu et al.l I2008T ) at z=0.062, the mean red- 
shift of the C I RS s ample, using the mass function of 
Uenkins et aD (p5ol . The WMAP5+SN+BAO con- 
straints combine the WMAP5 data with supernova data 
from recent survey s (lAstier et al.ll2006t 11 less 
IWood-Vasev et al.ll2007) and estimates o f baryon acous- 
tic oscillations from iPercival et alJ (|2007t ). The WMAP1 
and WMAP3 predictions straddle the CIRS results, 
while the WMAP5 and WMAP5+SN+BAO show ex- 
cellent agreement with the CIRS mass function. In fact, 
the mass function predicted by the WMAP5+SN-I-BAO 
parameters is almost indistinguishable from the best-fit 
CIRS mass function (green curve in Figure 1). 

Figure compares the CIRS cosmological constraints 
in the (fl m , 08 ) plane with constraints from WMAP5, 
WMAP5+SN-I-BAO, and cosmic shear. The green and 
yellow contours in Figure [2] show the 68% and 95% 
confidence levels for f2 m and as from the CIRS virial 
mass function. Blue, orange, and magenta contours 
show the 68% a nd 95% confidence levels from CFHTLS 
(|Fu et all 120081 ), WMAP5, and WMAP5+SN+BAO re- 
spectively. The 68% confidenc e levels of WMAP 3 and 
CIRS overlapped each other (Rin es et alJ [2007), but 
the best-fit value of WMAP3 suggested a smaller clus- 
ter abundance (Figure [T|). Both the WMAP5 and 
WMAP5+SN+BAO constraints have much smaller un- 
certainties than WMAP3. The CIRS constraints are in 
excellent agreement with these new, tighter constraints. 
Note that cosmic shear estimates of Q m and as have 
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Fig. 2. — Cosmological constraints from the CIRS virial mass 
function compared to other results. Green/yellow contours show 
68% and 95% confidence levels for f2 m and erg. Blue, orange, 
and magenta contours show the 68% and 95% confidence lev- 
els for CFHTLS, WMAP5, and WMAP5+SN+BAO resp ectively 
l|Fu et al.H200l IDunklev et~aHl200l IKomatsu e7al]|2008f ) . 

comparable statistical uncertainties with the CIRS con- 
straints. The cosmic shear constraints vary by approxi- 
mately the size of the 95% contours when either different 
measures of cosmic shear are used (e. g., two-point co rre- 
lation function versus aperture mass, IFu et al.ll200l or 
when applied to different samples (|Massev et al.l l 2007). 
These variations suggest that the uncertainties in cos- 
mic shear constraints are limited by systematic effects 
at present. The range of the cosmic shear constraints 
spans the range of the CIRS constraints, indicating broad 
agreement between these very different measures of large- 
scale structure. 

4. DISCUSSION 

We compare the recent revisions to cosmological pa- 
rameter estimates available from WMAP Five- Year data 
to the mass function of galaxy clusters in the local uni- 
verse. The differences between WMAP3 and WMAP 5 
have a large effect on the expected mass function. In 
particular, our determination of the cluster mass function 
found a larger cluster abundance th an expected for th e 
parameters estimated by WMAP 3 ijRines et all l2007f) . 
The revised parameters bring the expected mass func- 
tion into much better agreement with our data. 

Large velocity segregation between galaxies and the 
underlying dark matter could produce a significant off- 
set between the observed and true mass function (due 
to overestimating or underestimating cluster masses). 
Such large velocity segregation is not expected from 
state - of-the-art numer i cal simulations (iDiemand et all 
2004 iGao et all 12004 iFaltenbacher k. Diemandl 120061: 
Biviano et al.ll2006fl ; lEvrard et alJ (|2008D summarize the 
theoretical expectation as o" 9X1/ /cr£)jvf=1.0Q±0.05. With 
the WMAP5 parameters, modest velocity segregation 
(cgxy/ci3M~l-13^Q'Qg) is sufficient to produce agree- 
ment. When combined with constraints from supernovae 
and baryon acoustic oscillations, the implied velocity 
segregation is small {a gxy /auM ~1.05±0.05), consistent 
with expectations from simulations. If the WMAP5 and 
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CIRS results are correct, then cluster galaxies are indeed 
robust tracers of the velocity distribution (and dynamics) 
of the underlying dark matter. 

Similarly, the large excess specific energy required in 
the I CM to match WMA P3 results with X-ray observa- 
tions (|Evrard et al.ll2008l ) is also significantly reduced. In 
the parameterization 5' 8 =cr 8 (fi m /0.3) - 35 , WMAP3 in- 
dicated 58=0.69, while a larger value of 58=0.8 would 
imply that the ICM has comparable specific energy to 
the dark matter (perhaps 10% larger, which is within 
the range expected from simulations), while the im- 
plied ICM mass fractions would be consi stent with X- 
ray a nd Sunyaev-Zeldovich observations (|Evrard et al.l 
120081) . The WMAP5 results alone indicate 5 8 =0.75, 
while the WMAP5+SN+BAO results indicate 5 8 =0.80, 
exactly the value suggested above as a possible resolu- 
tion of the tension between WMAP3 and X-ray cluster 
observations. 



This agreement also supports the case for future large 
cluster survey s as p oten t ial co s mological prob e s(e.g. , 
1 l200l [Hul [20031: iVikhlinin et all [20031 : 
2004 iMantz et al.ll2008f ). We con- 



Haiman et al 
Maiumdar fc Mohrl 
elude that a large spectroscopic program to measure 
virial masses of an X-ray selected sample at moderate 
redshift would provide an independent measurement of 
the evolution of the mass function. If other cosmologi- 
cal probes provide tighter (and consistent) constraints on 
cosmological parameters than those achievable with clus- 
ter surveys, the data can be used to probe the dynamical 
history of clusters, e.g., by measuring the evolution of ve- 
locity segregation and the thermal history of the ICM. 



We thank Adrian Jenkins for providing his software 
for calculating the mass function. We thank Margaret 
Geller for suggesting that we write up this analysis. 



REFERENCES 



Adelman-McCarthy, J. et al. 2006, ApJS, 162, 38 
Astier, P. et al. 2006, A&A, 447, 31 
Bohringer, H. et al. 2000, ApJS, 129, 435 
— . 2001, A&A, 369, 826 

Bahcall, N. A. & Cen, R. 1993, ApJ, 407, L49 
Bahcall, N. A. et al. 2003, ApJ, 585, 182 
Becker, M. R. et al. 2007, ApJ, 669, 905 

Benatov, L., Rines, K. J., Natarajan, P., Nagai, D., & Kravtsov, 

A. 2006, MNRAS, 370, 427 
Biviano, A., Murante, G., Borgani, S., Diaferio, A., Dolag, K., & 

Girardi, M. 2006, A&A, 456, 23 
Diemand, J., Moore, B., & Stadel, J. 2004, MNRAS, 352, 535 
Dunkley, J. et al. 2008, ArXiv 0803.0586 

Ebeling, H., Edge, A. C, Allen, S. W., Crawford, C. S., Fabian, 
A. C, & Huchra, J. P. 2000, MNRAS, 318, 333 

Evrard, A. E. et al. 2008, ApJ, 672, 122 

Faltenbacher, A. & Diemand, J. 2006, MNRAS, 369, 1698 

Faltenbacher, A., Kravtsov, A. V., Nagai, D., & Gottlober, S. 2005, 
MNRAS, 358, 139 

Fu, L. et al. 2008, A&A, 479, 9 

Gao, L., De Lucia, G., White, S. D. M., & Jenkins, A. 2004, 

MNRAS, 352, LI 
Haiman, Z., Mohr, J. J., & Holder, G. P. 2001, ApJ, 553, 545 
Henry, J. P. & Arnaud, K. A. 1991, ApJ, 372, 410 
Hu, W. 2003, Phys. Rev. D, 67, 081304 

Jenkins, A., Frenk, C. S., White, S. D. M., Colberg, J. M., Cole, 
S., Evrard, A. E., Couchman, H. M. P., & Yoshida, N. 2001, 
MNRAS, 321, 372 

Kauffmann, C, Colberg, J. M., Diaferio, A., & White, S. D. M. 
1999a, MNRAS, 303, 188 

— . 1999b, MNRAS, 307, 529 



Komatsu, E. et al. 2008, ArXiv 0803.0547 
Loewenstein, M. 2000, ApJ, 532, 17 
Majumdar, S. & Mohr, J. J. 2004, ApJ, 613, 41 
Mantz, A., Allen, S. W., Ebeling, H., & Rapetti, D. 2008, ArXiv 
0709.4294 

Massey, R. et al. 2007, ApJS, 172, 239 

Nagai, D., Vikhlinin, A., & Kravtsov, A. V. 2007, ApJ, 655, 98 
Percival, W. J. et al. 2007, MNRAS, 381, 1053 
Pierpaoli, E., Scott, D., & White, M. 2001, MNRAS, 325, 77 
Rasia, E., Mazzotta, P., Borgani, S., Moscardini, L., Dolag, K., 

Tormen, G., Diaferio, A., & Murante, G. 2005, ApJ, 618, LI 
Reiprich, T. H. 2006, A&A, 453, L39 
Reiprich, T. H. & Bohringer, H. 2002, ApJ, 567, 716 
Riess, A. G. et al. 2007, ApJ, 659, 98 

Rines, K., Diaferio, A., & Natarajan, P. 2007, ApJ, 657, 183 

Rines, K. J. & Diaferio, A. 2006, AJ, 132, 1275 

Rozo, E. et al. 2008, astro-ph/0703571 

Schmidt, M. 1968, ApJ, 151, 393 

Spergel, D. et al. 2003, ApJS, 148, 175 

Spergel, D. N. et al. 2007, ApJS, 170, 377 

Stoughton, C. et al. 2002, AJ, 123, 485 

Tegmark, M. et al. 2006, Phys. Rev. D, 74, 123507 

Vikhlinin, A., Kravtsov, A., Forman, W., Jones, C, Markevitch, 

M., Murray, S. S., & Van Speybroeck, L. 2006, ApJ, 640, 691 
Vikhlinin, A. et al. 2003, ApJ, 590, 15 
White, M. 2002, ApJS, 143, 241 
Wood-Vasey, W. M. et al. 2007, ApJ, 666, 694 
Zwicky, F. 1933, Helv. Phys. Acta, 6, 110 
— . 1937, ApJ, 86, 217 



